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We address the forces exerted by the electromagnetic field emitted by a planar fluctuating source 
on dielectric particles that have arose much interest because of their recently shown magnetodielec- 
tric behavior. In this context, we analyze as a particular case the modification of the Casimir and 
Van der Waals forces. We study the effect of the source coherence length as well as the interplay 
between the force from the radiated field and that from the electric and magnetic dipoles induced 
on the particle. This allows a control of these interactions as well as of the weight and interference 
effects between the fields from both kinds of induced dipoles, in particular when large changes in 
their differential scattering cross section occur due to Kerker minimum forward or zero backward 
conditions; thus opening new paths to nanoparticle ensembling and manipulation. The influence of 
surface waves of the source is also studied. 
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Magnetodielectric particles have recently attracted 
much attention due to their exotic properties as scatter- 
ers and nanoantennas [1-9]. It has been shown [10, 11] 
that some dielectric particles behave in this way exhibit- 
ing coupled electric and magnetic dipoles induced by the 
illuminating light field. It is then of great interest to 
study their response to stochastic radiation forces, specif- 
ically those due the field of a planar fluctuating source 
[12]. Many previous studies have dealt with this subject 
for atoms or non-magnetic nanoparticles with its appli- 
cation to delta-correlated thermal sources and blackbod- 
ics in connection with the Van der Waals (VdW) and 
Casimir-Polder (C-P) interactions [13-15]. 

In this work we deal with a more general kind of sta- 
tistical sources, namely those that are spatially partially 
coherent and of the wide variety of those statistically ho- 
mogeneous and isotropic [12, 16]. Their emission excites 
electric and magnetic dipoles of the particle in its near 
field that may be considered as a secondary source whose 
radiation interacts with the primary fluctuating source; 
this giving rise to a new total force resulting from both 
the action of the field radiated by the primary source and 
that from the field emitted by these dipoles. At thermal 
wavelengths, the interaction from this secondary source 
constituted by the particle induced dipoles is interpreted 
as a Liftshitz force [17], which at zero temperature be- 
comes either those derived by VdW and C-P [18], de- 
pending on the distance and use, or not, of quasistatic 
formulations. However, in our study, first the optical 
wavelengths are such that Hui/kT » 1 and hence Planck 
energy becomes similar to that of the vacuum fluctua- 
tions: hu)[\ + 1 / (exp(huj / kT) — 1)] s=s i/kj; and second, 
due to the magnetic response of the nanoparticle, more 
forces come into play in addition to those that keep an 
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analogy with these above quoted, thus allowing a larger 
number of degrees of freedom of relevance for particle 
ensembling and manipulation. We do not restrict here 
to thermal wavelengthts, but rather consider a range of 
infrared and optical frequencies whose choice depends on 
the particle size. Then a particular case of our study is 
that analogous to those dipole forces induced by sponta- 
neous electromagnetic field fluctuations. While the vac- 
uum forces can become relatively negligible by conve- 
niently manipulating the power intensity of the source, 
we shall show that in other possible configurations and 
experimental designs they may predominate. 

The geometry considered in this letter consists of two 
half-spaces. The lower one (z < 0) is occupied by the 
source with its polarization currents and it will be de- 
noted as 1; whereas the upper one (z > 0), denoted as 2, 
is free space and contains the particle. 

The Cartesian components of the force on a magne- 
todielectric dipolar particle is the sum of an electric, 
magnetic and electric-magnetic dipole interference parts 
which are expressed in terms of the first electric and mag- 
netic Mie coefficients a\ and b\ as [19] 

F i (v)=F?(r) + Fr(r)+Fr m (r) 
= S -^Re{(a e E*(r)d i E j (v))} 

+ ^lRe{(a e H*(r)d i H j (r))} 
7k 4 

-e e 2j ^Re{(a* e a m )(E*xK) i }, (1) 

where i, j = 1,2,3, ei = e\ + ie" and \ii = /ij + i/i" 
(/ = 1,2) are the permittivity and susceptibility of the 
medium embedding the particle, respectively, in our case 
being vacuum; and Z = d ^o/xi/eoei- a e and a m stand- 
ing for the electric and magnetic polarizability of the par- 
ticle, respectively. a e = ifp-ai, a m = i^j^h; Ei(r) is 
the total electric vector at frequency uj at any point of 
the half-space z > 0, hence at the position of the particle, 
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i.e. at r = r , it will be 

Ei(T ) = E™ c (r ) + Ef(r Q ) + £f (r ) 

= ET C + noW^jPj + inoWjGKmj . (2) 

In the second line of this equation we have omitted the 
explicit dependency on space and frequency for brevity. 
The associated magnetic field can be obtained directly 
from Maxwell's equations. In Eq. (2) the last two terms 
are the electric fields emitted by the particle induced 
dipolcs (pi = eoot e El nc and rrii = a m Hl nc ) after mul- 
tiple reflections at the plane z = [20] and hence con- 
nect the constitutive properties of the source and particle 
through the reflection Fresnel coefficients r sp at z = 
and the polarizabilities, being described by G^j" 1 . On 
the other hand, the first term E\ nc represents the electric 
field incident on the particle after being emitted from the 
primary source placed at z < 0, and is defined through 
the Green's function Gf, which includes the transmis- 
sion Fresnel coefficients t StP from z < into z > 0. 
Gfj may be written as a superposition of plane waves 
with wavevector k^Si = fc (sj_, s z i ) with Sj_ = (s x ,s y ), 
s 2 z i = EifM — Sj_,(i = 1,2), and fco = 2ir/\ — cu/c. Thus, 
in terms of the polarization currents one has 

Er c (r ) - W>^ 2 / Gfj {t , r', co) P 3 (r', u ) dV. (3) 

In vacuum (ff 2 = A*2 = 1), G p p decays exponentially 
with the distance in the evanescent wave region (s± > 
1) and is oscillatory in the radiative one (s± < 1) [21]. 
For a single dipole, this electric field is calculated from 
Eq. (3) using Gfj and PiS(r' — r ) instead of Gf 3 - and 
Pj(r'). In a similar way, the electric field E m generated 
by the magnetic dipole can be calculated from the Green 
function associated to the magnetic field of the electric 
dipole H p and making the interchange r s <R- r p in the 
reflection Fresnel coefficients. Some details about these 
Green's functions can be found in e.g. [22, 23]. Note 
that for p-polarization, both t p and r p support surface 
plasmon polaritons (SPPs) when e' 2 < 1. 

Now we define the cross-spectral density tensor 
of the source polarization as (ri,r 2 ,w) = 

(P* {v\,Lj)Pj (r 2 ,w)). Wc shall address the wide vari- 
ety of non-local statistically homogeneous and isotropic 
sources [12]) for which 

W<>P (n,r 2 ,a;) = S^\u)t^\\v x r 2 |, W ). (4) 

S p (ui) denotes the power spectrum of the source and 
Hij(ri,r2,u)) is the spectral degree of coherence [12]. A 
special particular case of these sources are those ther- 
mal and blackbodies widely studied. On inserting Eq. 
(2) into (1) and taking the statistical homogeneity of the 
source into account, one obtains the total force on the 
particle. For these sources only the force along the z- 
axis is different from zero. We consider mutual incoher- 
ence between the particle electric and magnetic induced 



dipoles, i.e., (p*m 3 ) — [24]. F\ and Fi will denote the 
total forces due to the above mentioned contributions 
of the primary fluctuating source in z < and to the 
secondary source field from the particle induced dipolcs, 
respectively. 

We study the effect of the magnetodielectric properties 
of the particle in the near infrared. This is a semicon- 
ductor sphere; its anomalous scattering properties have 
recently received a great deal of attention, both theoreti- 
cally and experimentally [3, 11, 25, 26]. In particular, for 
each incident plane wave component its scattered inten- 
sity in the backscattering direction is zero, [first Kerker 
condition (Kl)] when Rea e = Rea m is fulfilled. Also for 
each of such plane wave components impinging the parti- 
cle, the forwardly scattered intensity becomes close to a 
non-zero minimum [second Kerker condition (K2)] when 
Rea e = — Rea m . In both cases: Ima e = Ima m [25]. 

Thus, we address a Si sphere of radius a = 230nm, 
the spectrum of the incident light being in the range of 
1.2 — 2/jm. At these frequencies the total cross section 
of the particle is fully determined by the Mie coefficients 
ai and bi [10], this justifies the use of Eq. (1). 

We assume a Gaussian degree of co- 
herence of the source, therefore the cor- 
relation function reads (ri,r2,u;) = 

5( p ) Mcxp (- (|n - r 2 |) 2 /2a 2 ) ^•/(2 7 r) 3 / 2 C 7 3 , 

being the coherence length of the source and 
£( P V) = 5( p )(w)/(27r) 3 / 2 cr 3 representing the 
normalized spectrum. Notice that for a — > 0, 
W$p{r U T 2 ,w) = S^MJdn - r 2 |, W )% then 
the source becomes 8-correlated, like in e.g. a thermal 
one following the fluctuation-dissipation theorem [24]. 
The source considered here will have an Au interface 
at z = 0, hence supporting surface plasmons polaritons 
(SPPs) in the spectral range under consideration, 
and thus enhancing the near field forces [27]. After 
performing all the integrations of the form of Eq. (3) 
for the primary source and for the induced dipoles, (a 
long but straightforward work, some of whose details 
are shown in the supplementary information), one sees 
that (£f*(r )£r(ro)) = (Hf*(r )Ef(r )) - 0. This 
will be relevant when we discuss the particle Kerker 
conditions. Fig. 1 shows the logarithm of the force. 
This representation aims to clarifying its drastic changes 
of sign. All the results of this paper will be normalized 
to the spectrum of the source in order to see the relative 
weight of each force component. The first horizontal row 
shows the force F± due to the field impinging from the 
primary statistical source at z — 0. The second horizon- 
tal row represents the force Fi from the secondary source 
constituted by the electric and magnetic dipoles induced 
on the particle. The fluctuating source coherence length 
a is first assumed to be zero. The inset shows the be- 
havior of the polarizability for the range of wavelengths 
considered. This helps to understand the color plots. In 
Figs. 1(a) and 1(b) we see a line separating the gradient 
and scattering forces. For a statistically homogeneous 



3 



Log (Fe) (arb. units) 



Log (Fm) (arb. units) 



Log (Fern) (arb. units) 




FIG. 1. (Color online). From left to right: Normalized elec- 
tric (F e ), magnetic (F m ) and interference force (F em ) from a 
5-correlated source. The first horizontal row shows the part of 
the force due to the field from the fluctuating primary source 
at 2 = 0, [first term in Eq. (2)]. The insets exhibit the behav- 
ior of the electric, magnetic and electric-magnetic interaction 
force components, respectively, normalized to a 3 . The second 
horizontal row represents the force from the secondary source 
constituted by the particle induced dipoles, [second and third 
terms in Eq. (2)] 



source, the gradient force (proportional to Rea e (E*Ei)) 
is governed uniquely by the evanescent modes and is 
negative for a particle with Rea e > [27], hence, it 
exponentially decays with the distance z to the source. 
On the other hand, the scattering force (proportional to 
lma e Im(Ej diEj}) is positive, i.e. pushing, and constant 
for any ro. As the wavelength grows, Rea e > Ima e , 
[see the inset in Fig. 1(a)], the gradient force dominates 
even at distances larger than A, where the evanescent 
modes do not contribute. This is a remarkable new 
feature of this kind of particles. Fig. 1 (c) represents the 
force component F^~ m due to interference between the 
particle induced electric and magnetic dipoles. In the 
near field this is almost repulsive for any wavelength, 
however, at distances larger than the wavelength, where 
the Poynting vector S = Re(E* x H)/2 is independent of 
the distance, we have a zone where this force is negative 
(the arrow indicates that zone). This kind of action 
is known as a pulling force [28-30], and its interest 
has increased in the last years. This last plot shows 
the relevant role of the magnetodielectric behavior of 
these particles in this respect, although in this latter 
specific case when the two other components: electric 
and magnetic, are added this pulling effect becomes very 
small, (the electric-magnetic dipole interference force at 
A ~ 1.47/zm and z > 2.5A is one order of magnitude 
less than either Ff or F™), however by manipulating 
the fluctuating source, a suitable power of the emitted 
electromagnetic field is obtained [15] so that a tractor 
light field appears, even in the far zone. 

Concerning the force i*2 fom the particle induced 



dipoles, we observe in the second row of Fig. 1 that 
this force exponentially decays with the distance z to the 
primary source plane z = 0, and its sign depends on that 
of the particle polarizability; nevertheless, the oscillating 
behavior of the Green function due to propagating plane 
wave components manifests in this force. We also observe 
that it is six orders of magnitude larger than its counter- 
part F\ from the primary fluctuating source, at least at 
subwavelength distances z. We shall later discuss this. 

To get a deeper understanding, Fig. 2 represents F 2 
for some selected wavelengths and for two different source 
coherence lengths: a — and a = A/4. For an statisti- 
cally homogeneous source, the relationship between the 
electric and magnetic cross-spectral density tensors is [31] 

e <^(r 1 ,«)£, J (r a ))= A io(fl?(ri,a;)flJ(r 2 )) ) (5) 

hence, in the near field, when the first or the second 
Kerker condition holds, one has Ff = F™ and Ff = 
— Ff 1 , however, in the far zone Ff = FJ™ for any value of 
ro. For the Si particle addressed, the Kerker conditions 
are fullfilled at Ai ~ 1.825/im and A2 — 1.53/im, (see 
Figs. 2(a) and 2(b) and [25]). We can also see in the inset 
of Fig. 1(b), that in the range of A = 1.6 — 1.65/xm there 
is a peak in the imaginary part of a m which predominates 
over all other a parts. The black-dashed-dot line in Fig. 
2 represents the force in this peak. In that case, the 
magnetic force F£\ is one order of magnitude larger than 
Fgi and Fff, hence the total force on the particle is 
governed by this magnetic force. This effect, due to the 
dielectric particle magnetic response to the light field, 
constitutes one of the main results of this paper. 




-10 
10 

8 

«T 6 

1 1 

•S 2 

3. 

I- 2 
o-« 

_l -6 
-8 
-10 



<b) 



o=0 



***** 



=1.825 (k^ 
=1-53 (k ) 



■ A=1 .65 

■ 1=1 .25 



(<J) 


o=W4 















(e) 


CT=V4 









(f) 



).5 1 1.5 2 2.5 3 0.5 1 1.5 2 2.5 3 0.5 1 1.5 2 

zA zlX zlX 



FIG. 2. (Color online). Normalized Ff, and Ff m versus 
the distance from the plane of the source (in wavelength units) 
for some values of the wavelength (in (im). The two Kerker 
conditions occur at Ai, (Kl), and A2, (K2), respectively. 

We now address the influence of the coherence length of 
the source. This will establish the differences between the 
mechanical action of partially correlated sources and that 
from e.g. thermal sources and blackbodies. The spec- 
tral degree of coherence in k— space is exp[— (kas±) 2 /2], 
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hence, it acts a low-pass filter being maximum for a = 
(i.e. when the source is ^-correlated) . Because of this 
fact, the evanescent modes present two such filters: the 
first is due to the own nature of these evanescent modes 
while the second stems from the spatial coherence of 
the source. The shape of Figs. 2(d)-(f) is similar to 
that of Figs. 2(a)-(c), shifted by a distance Az ~ 0.5A, 
therefore for a > A the force is solely due to the non- 
conservative (scattering) force and to the interference 
force Ff m , which becomes constant and positive or neg- 
ative depending on the wavelength. It is worth pointing 
out that the price paid on increasing the coherence length 
is expensive, because at the same time there is a reduc- 
tion of the force strength by various orders of magnitude 
[cf. e.g. the forces shown in Figs. 2 (a) and 2 (d)]. 
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FIG. 3. Color online). Normalized F 2 e , F 2 m and F 2 em versus 
distance 2 (in units of wavelength) from the exit plane 2 = 
of the source for different values of the wavelength (in /jm). 
The two Kerker conditions are fulfilled at Ai, (Kl), and A2, 
(K2), respectively. 



near-field z < A is much larger than in Fig. 2, thus, the 
effect of the mechanical action F2 of the field emitted 
by the particle induced dipoles substantially dominates 
over that Fx of the field that is due solely to the stochas- 
tic source. Nevertheless, as the distance z grows, all the 
fastly oscillating components, electric, magnetic and that 
of interference, of this force F 2 rapidly tend to zero, and 
hence is the force F% from the primary source the one 
that dominates. As follows from the calculation of E p,m 
and H p,m , the cross spectral density tensor of the elec- 
tric and magnetic dipoles now are not equal; therefore, 
and although at first sight it could seem that similar re- 
lationships between the electric and magnetic forces, in 
Kerker conditions, are fulfilled like for F\ from the pri- 
mary source, in fact they are not. 

The role of the coherence length in this case is exactly 
the same as in Fig. 1; the magnitude of the force de- 
creases as a grows. Future work should find a minimum 
value of a for which the Casimir-Polder force predomi- 
nates over the contributions discussed here. 

In summary, we have shown, that the mechanical ac- 
tion on a magnetodielectric small particle from a par- 
tially coherent fluctuating source exhibits important new 
effects at distances shorter than the wavelength. In par- 
ticular, the magnetic induced dipole and its interference 
with the electric dipole create a landscape of forces com- 
pletely different to that previously studied in connection 
with Van der Waals, Casimir, Liftshitz and rest of ra- 
diation forces both in and out of thermodynamic equi- 
librium. In this respect, Kerker conditions, as a result 
of the exotic scattering properties of these particles, in- 
troduce new relationships in the balance of these new 
forces over the traditional purely electric forces. Further 
experiments should be stimulated by these new effects, 
including consequences of Fano resonances. 
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